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Transverse-momentum spectra in hard processes are typically described either in terms of intrinsic 
transverse momentum of partons, or in terms of perturbative radiation. The relation between these 
descriptions is discussed for the example of semi-inclusive deep inelastic scattering, with special 
focus on the angular distribution of the observed hadron. This involves nontrivial theoretical 
issues, such as the proper definition of transverse-momentum dependent parton distributions, and 
has practical consequences for the description of pj spectra in phenomenology. 



X 
J3 



LIGHT CONE 2008 Relativistic Nuclear and Particle Physics 
July 7-11, 2008 
Mulhouse, France 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Transverse momentum of partons: from low to high pT 



Markus Diehl 



1. Motivation 

This talk is concerned with the transverse-momentum dependence of particles produced in 
scattering processes involving a hard momentum scale. Specifically, I will discuss semi-inclusive 
deep inelastic scattering, ep e + h + X, where the momentum of the hadron h is measured. 
Using crossing symmetry, the results can be carried over to Drell-Yan or W or Z production in pp 
and pp collisions, as well as to hadron pair production in electron-positron annihilation, e^e^ —>■ 
/ji +/i2 +X. Details can be found in the recent paper [jl|]. 

The transverse-momentum spectrum of produced particles may be regarded as a basic feature 
of the final state. Even in the simple processes just mentioned, its investigation reveals a number 
of nontrivial aspects of QCD dynamics. There are two theoretical frameworks to describe the dis- 
tribution of a suitably defined transverse momentum qj in the final state. The description sketched 
in Fig. [T]a is based on the "intrinsic transverse momentum" of partons within a hadron and uses 
transverse-momentum dependent (i.e. unintegrated) parton densities and fragmentation functions. 
This description can be used for qj <C Q, where Q is the virtuality of the photon or electroweak 
boson, which I assume to be large throughout this talk. The description represented in Fig.[I]b uses 
collinear (i.e. kj integrated) parton densities and fragmentation functions, generating finite qj by 
perturbative radiation of partons into the final state. This description is adequate for qj ^ M, 
where M stands for a generic nonperturbative scale. In the following I refer to the two descriptions 
as "low-qr" and "high-^r", respectively. It has long been known that both mechanisms give rise to 
a nonzero cross section for longitudinal photon polarization and to a nontrivial dependence on the 
azimuth of 




Figure 1: The low-^j- description (a) and the high-g'T- description (b) for the transverse-momentum distri- 
bution of the produced particle h in semi-inclusive deep inelastic scattering, ep^e + h+ X. 



2. Insight from power counting 

It is natural to ask how these two descriptions are related to each other. A first answer can 
be obtained from a careful look at the power counting in the region of intermediate transverse 
momentum, M <^ qj ^ Q, where both approaches can be applied. The low-qj approach starts with 
an expansion in the small parameter qr/Q and involves coefficients depending on M/qj, which for 
qr ^ M can be further expanded in M/qj. For an observable F with mass dimension —2 one thus 
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By contrast, the high-^7 approach first expands an observable in M/^y, with coefficients that for 
intermediate qj can be further expanded in qj/Q: 
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The simultaneous validity of both approaches in the region M <^qj <^Q implies Z„ ,t = h^^n- The 
first index « > 2 in each expansion characterizes the twist of the corresponding calculation. In 
practice only terms with n = 2 and possibly n = 3 can actually be calculated. 

For observables with nonzero h;! = /22,2> the leading terms in the two calculations coincide 
for intermediate qT, where they provide complementary descriptions of the same physics. One 
may then try to construct a smooth interpolation between the two descriptions that is valid at all qj. 
There are, however, observables with I22 = ^2,2 = 0, where the leading term /2,4 of the low-^r result 
is distinct from the leading term /i2,4 of the high-g^ result. With both calculations only performed at 
leading-twist accuracy, one can then add their results at intermediate qT without double counting; 
which of them is more important at given qj depends on the relative size of qj /Q and Mjqj. We 
will encounter examples for both cases in section H. 



3. Structure functions for semi-inclusive deep inelastic scattering 

To describe the kinematics we use the standard scaling variables x and z, the inelasticity y, the 
photon virtuality Q, the scaled transverse momentum qj = Ph±/z of the produced hadron in the 
Y*p cm., and the azimuthal angle <p between the lepton and hadron planes in that frame. Precise 
definitions are given in The unpolarized cross section can then be parameterized in the form 

do{ep ehX) 



d<^ dqfdxdydz 
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Ft + sFt + V2e(l + e) cos ^ + £ cos 20 F'°' ^0 



(3.1) 



where the ratio of longitudinal and transverse photon flux is given by e = (1 —y)/{l —y + y'^/2) 
in the Bjorken limit. The semi-inclusive structure functions F depend on x, z, q^, Q^, and the 
subscripts T and L are respectively associated with transverse and longitudinal photon polarization. 



3.1 High-^r calculation 

The high-^r calculation gives the structure functions as convolutions 
1 ^ r^dx 
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in longitudinal momentum fractions, where /J and Dj respectively denote the usual unpolarized 
collinear distribution and fragmentation functions, and the K'-l are perturbatively calculable hard- 
scattering kernels. Expanding these kernels in qj / Q, one readily obtains the result of the high- 



qj mechanism for intermediate qr- At order a,, all structure functions in (3.1) contain a term 
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proportional to /i (x) D\{z) log (Q/qr) after this expansion. Higher orders provide terms going like 
a" log ^'"^^(2/^7-) with m<n.To obtain a stable perturbative result in the region where log(2/^7-) 
is large, these logarithms should be resummed to all orders. We will come back to this in sect. 

3.2 Low qr 

In the low-^77- description, factorization is fairly weU understood at twist-two accuracy ^ ^ 
and leads to a representation of the form 



d^Pj d% d% 8{pj-kj + lj + qj)w,„{PTM)f PT)D'{z,k^)U{l^) (3.3) 



with known functions w... , where for simplicity I have omitted a hard factor representing 01^ correc- 
tions from virtual graphs. In addition to transverse-momentum dependent distribution and fragmen- 
tation functions f'{x,p^) andZ)'(z,^y), the expression ( |33| ) contains a soft factor U{1^) describing 
soft gluon exchange between partons moving in the direction of the target and partons moving in 
the direction of the observed hadron h. At twist-three accuracy, soft gluon exchange has not been 
analyzed, so that we do not have a full understanding of factorization for structure functions go- 
ing like l/Q. However, there are detailed calculations at tree level [M, 0], which give results very 



similar in form to ( |3.3[ ) without the factor U{lj). 



To evaluate ( |3.3[ ) for intermediate qj, one notes that for qr ^ M at least one of the momenta 
Pt, kj or It has to be large. In this region, the corresponding factor can be calculated using collinear 
factorization, with the large transverse momentum generated by perturbative parton radiation as 
illustrated in Fig. ^ As shown in general symmetry considerations allow one to determine the 
power behavior at high pt for the different parton densities parameterizing the spin and momentum 
dependence of the quark distribution in a proton. In particular one finds 



hf{x,pl) 



'=i,q-. 



Pt i=q,q,g Pt 

— Us 2^ [collinear twist-three distributions <E> K'^] , 
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(3.4) 



where K\, K^', and A'j are calculable hard-scattering kernels and (g) denotes a convolution in lon- 
gitudinal momentum fractions. The Boer-Mulders function hf^ describes the distribution of trans- 





(a) 



(b) 



Figure 2: Example graphs for the calculation of an unintegrated parton distribution at high pr in terms 
of a collinear distribution f{x) and a hard-scattering subprocess. The double line represents an eikonal 
propagator, which comes from the Wilson line Pe\p[—ig J^dXn ■A{^ +A«)] in the definition of theparton 
density 
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versely polarized quarks in an unpolarized proton, whereas is a distribution of twist three (in- 
volving one good and one bad light-cone component of the quark field). Explicit calculation of 
/i {x,Py) and f^{x,pj) reveals the rapidity divergences discussed in They have the form 

111 pi 

with {l-p)- = —t^— (3.5) 



{l-p)-n n- (i_p)+^!^{^i_p)- 2{l-p) 

and come from the Wilson line in Fig. or equivalently from the gluon propagator in Fig. ^ 
if one uses the gauge n ■ A = where the Wilson line is unity. In a parton density for a fast 
right-moving proton, the loop variable Z+ is integrated down to its lower kinematic limit To 
avoid a logarithmic divergence in ( |3.5| ) one must hence keep nonzero (complications arising for 
spacelike n are discussed in ^). The light-cone gauge A+ = or a purely lightlike Wilson line 
cannot be used in this context. It should be instructive to investigate how this affects formulations of 
QCD based on light-cone gauge. Note that in the context of light-cone quantization, configurations 
with (/ — ;?)+ = in Fig. ^ correspond to zero modes of the gluon field. 

Keeping finite in (3^) cuts off the region of negative gluon rapidities, where (/ — p)+ — > 



and {I — p)^ — > oo. This is physically reasonable, since fast left-moving gluon modes should not be 
included in the parton distribution of a fast right-moving hadron. The dependence of unintegrated 
parton distributions on jnT is described by the Collins-Soper equation [|5]], whose solution can 
be written as the product of an jn^ independent initial condition and a Sudakov factor. 

Power laws analogous to those in are obtained for the fragmentation functions Di(z,^y), 
D^(z,^y), and the Collins function H^{z,kj). Together with the perturbative expression for U {Ij) 
at high Ij one can then determine the behavior of the semi-inclusive structure functions for M -C 
qT *C Q. The terms going with jnT in the distribution and fragmentation functions give rise 
to a log(2/^7-) in the structure functions, which we also encountered in the high-^r calculation. 
The Collins-Soper equation allows one to resum such logarithms to all orders and is at the origin 



of the CSS formalism [10], which plays a prominent role in collider phenomenology. The need to 
keep jnT finite in (3.5) is thus not a mere technicality but has practical imphcations for physical 
observables. 



4. Comparing the low- and high-^^ calculations 

The behavior for M ^ ^7- ^ 2 of the unpolarized structure functions obtained in the low- 
and high-^r calculations is given in Table [ij Explicit evaluation shows that for Fj the results 
of the two calculations exactly coincide; this is an example of the case where at intermediate qj 



the leading term in the expansions ( [2. 1[ ) and ( p^ ) goes with I22 = /J2,2- The agreement between 
the two calculations can be seen at the level of diagrams: roughly speaking, the graph of Fig. |I|b 
corresponds to the graph in Fig. ^a if the gluon moves fast in the direction of the hadron h, and to 
the graph in Fig. ^ if the gluon moves fast in the direction of the target. 

The leading term in the low-^y result for f'=°s20 involves the Boer-Mulders and Collins func- 
tions. This structure function provides an example for the case where I2.2 = /J2,2 = and where 
the leading contributions /2,4 and /j4 2 in the two calculations are different and can be added at 
intermediate qj. From power counting it is clear that the 1/2^ behavior obtained in the high-<77' 
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Table 1: Behavior of semi-inclusive structure functions in the intermediate region M <^ qr ^ Q- 



calculation for both P'^°^'^^ and Fi corresponds to twist-four contributions in the low-^y frame- 
work, whose computation is well beyond the state of the art. As a consequence, one cannot invoke 



the CSS method [ |10| ] to resum the logarithms in Q/qj that appear in the high-^r result. Terms 
going hke 1/2^ are obtained if one calculates Fl and F^°^^^ in the parton model [Q], considering 
the graph in Fig. |I]a with only the functions f\{x,p\) and Di(z,/c^). However, the results do not 
match with the high-<77' calculation at intermediate qj and can hence only be regarded as partial 
evaluations of the complete (unknown) twist-four terms. 

If we perform the twist-three calculation for F'^°^ ^ at low qj using the tree-level result of 
supplemented with the soft factor U{lj) of the twist-two factorization formula (^), we obtain 
agreement with the high-^r result for intermediate qj except for a missing term proportional to 
/i(x)Di(z). Such a term comes from kinematics where both the plus- and the minus-momentum 
of the gluon in Fig. ^ is negligible. This shows that, if a proper factorization formula for twist three 
can be established, the soft-gluon sector will have to be treated with particular care. 




D(z,kT) 




f(x,PT) 



Figure 3: Example graphs for the low-qr calculation in the region qj ^ M, where a factorized description 
as in Fig. ||is valid for the fragmentation function (a) or the parton distribution (b). 



5. Summary 

The descriptions of transverse-momentum spectra based either on intrinsic transverse mo- 
mentum of partons or on perturbative radiation are not disconnected. They can be related in an 
intermediate region M <^ qj <^ Qby describing transverse-momentum dependent distribution and 
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fragmentation functions themselves in terms of perturbative radiation, as indicated in Fig. 0. Under- 
standing the connection between the two approaches for a given observable enables one to devise 
descriptions that may be used in the full region of qr- As shown in |jl|], the interplay of the two 
mechanisms has also nontrivial consequences in observables that are integrated over qj. 

A varied picture arises already for the angular distribution of the measured hadron in unpo- 
larized semi-inclusive scattering, with the results obtained at leading-power accuracy in the low- 
and high-^r calculations coinciding for some observables but not for others. An even richer phe- 
nomenology emerges if one includes polarization effects [[T], [TT]]. 

The calculation of unintegrated parton densities or fragmentation functions at perturbatively 
large transverse momenta leads to divergences from gluonic zero modes if naively performed in 
light-cone gauge. Proper regularization of these divergences physically ensures that left-moving 
gluon modes are not included in distribution or fragmentation functions for right-moving hadrons 
and provides a powerful method for resumming logarithms of Q/qr into a Sudakov factor. It 
remains to be understood how this physics can be treated within light-cone quantization. 
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